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Quasiparticle Heat Transport in Mixed State of High-T c Superconductors 
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The field dependence of low-temperature thermal conductivity k(H) observed on cuprates is 
explained by calculating n(H) microscopically. The heat current carried by low-lying quasipar- 
ticles around a vortex core decreases with H due to its depletion by the if-induced moment 
centered at a core for an underdoped case. This leads to a common and consistent picture on 
an "anomalous core". k(H,T) is found to be a useful tool for probing quasiparticle structures 
in the mixed state in general. 

KEYWORDS: thermal conductivity, vortex lattice state, antiferromagnetism, Bogoliubov-de Gennes theory, 
extended Hubbard model 



Much attention has been focused on high-temperature 
cuprate superconductors since their discovery over 15 
years ago. Although no consensus has been reached yet 
as to the pairing mechanism to induce high T c , we are 
gaining deep basic knowledge regarding unconventional 
d-wave superconductivity in general. Recently, a series of 
experimental and theoretical studies have revealed sig- 
nificant phenomena associated with a quantized vortex 
in high-T c superconductors. A mere few Tesla field has 
discovered hidden superconducting ground state prop- 
erties that could not be observed by conventional bulk 
measurements. Namely, a magnetic field applied to the 
c-axis induces magnetic moments exclusively around a 
vortex core. This is evidenced by neutron scattering ex- 
periments 1-4 which show an enhancement of the static 
incommensurately modulated moments under a field in 
La2-a;Sr a; Cu04+ ( 5 (LSCO). A muon spin resonance ex- 
periment 5 also shows an enhanced magnetic activity un- 
der a field. The antiferromagnetism (AFM)-induced vor- 
tex state is reproduced by some theoretical works. 6-10 

In addition to these experiments which probe the mo- 
ment increment averaged over an entire sample, other 
local probes reveal that magnetic moments exclusively 
appear locally around a vortex core. A scanning tunnel- 
ing microscopy (STM) experiment 11 on Bi2212 directly 
reveals a checkerboard pattern of the local density of 
states (LDOS) with a four-atomic-site period around a 
core, which strongly matches the neutron results with 
eight-site spin modulation. Site-selective NMR experi- 
ments on YBCO 12 and Tl 2 Ba 2 Cu0 6+2; 13 ' 14 also demon- 
strate that the relaxation time at the core site becomes 
longer than those at other sites away from the core, im- 
plying that LDOS at the core site decreases with H. 
These experiments collectively point to the "anomalous 
core" in cuprates. This picture allows us to understand a 
long standing puzzle ( "empty core" ) of the absence of the 
zero-energy DOS (ZEDOS) peak expected for a d-wave 
pairing vortex as observed by STM. 15-17 

Recent thermal conductivity measurements 18 ' 19 on 
LSCO under H are interesting in this respect because 
heat current at low T is carried by quasiparticles associ- 
ated with vortices. Thus, thermal conductivity provides 
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yet another important information concerning LDOS 
around a vortex. Two independent studies by Sun et 
al. 18 and Hawthorn et al. 19 yield essentially the same 
result: According to Sun et al. (1) for an overdoped re- 
gion (x = 0.17 and x = 0.22), thermal conductivity 
k(H)/k(Q) increases as a function of H at low T while 
it decreases as T increases. At an intermediate tempera- 
ture T*, k(H)/k(Q) becomes almost constant as H in- 
creases (Fig. 2 in ref. 18). (2) For a sample with an 
intermediate doping x = 0.14, k(H)/k(Q) decreases at 
low T. This decreasing tendency is strengthened as T in- 
creases (Fig. 3(c) in ref. 18). (3) For underdopings with 
x = 0.10 and x = 0.08, k(H)/k(Q) under a fixed field 
is minimum at a finite temperature (see Figs. 3(a) and 
3(b) in ref. 18). (We define this characteristic tempera- 
ture as T min .) The essential features of these results are 
also reproduced by Hawthorn et al.; 19 in particular, at 
the lowest T(= 60 mK), k(H)/k(0) decreases (increases) 
for underdoped (overdoped) samples (Fig. 2 in ref. 19). 
They attribute this behavior to "field-induced thermal 
metal-to-insulator transition" . Similar behaviors are also 
observed in YBCO. 20 

Here, we are going to explain these seemingly complex, 
yet rich behaviors by calculating k(H, T) microscopically 
in the mixed state for d-wave superconductors. 21 We 
assume the possibility that underdoped cuprates have 
an intrinsic tendency towards AFM instability. This as- 
sumption is completely justifiable both experimentally as 
mentioned above and physically because it turns out that 
the induced moment appears exclusively at the core sites 
where superconductivity is locally weakened, revealing a 
"true" ground state when removing superconductivity. 10 

Before going into detailed computations, we provide 
a clear view for understanding these phenomena. In the 
overdoped region away from magnetism, a standard d- 
wave vortex picture must be valid. n{H) is basically an 
increasing function of H at low T because of the growth 
of ZEDOS. At high T more than T* , n(H) decreases be- 
cause of the inability of the ZEDOS at the core to con- 
tribute to k in high T whose region increases with H (see 
the detailed exposition in ref. 22). Towards underdop- 
ings, the system exhibits AFM instability exclusively at 
the vortex core, leading to the removal of ZEDOS. Thus, 
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n(H) decreases with H . This instability is further ampli- 
fied as H grows since the AFM moment is enhanced. 

We start with an extended Hubbard model on a two- 
dimensional square lattice, and introduce the mean fields 
ni,<7 = ( a l a a i,a) at the i-site, where o is the spin index 
and i = (i x ,i y ) and Ae,i,<? = { a i,-aO-i+e,a) ■ We assume 
a pairing interaction V between nearest-neighbor (NN) 
sites. This type of pairing interaction gives ci-wave su- 
perconductivity. Thus, the mean-field Hamiltonian under 
# 10,23 is given by 

ft = - XI ti>j a l<r a i>° + U ^2 n i-^ a U a i,T 

i,j,<7 i,cr 

(1) 

where a\ a (oi, CT ) is the creation (annihilation) operator, 
and i + e represents the NN site (e = ±x, ±y). The trans- 
fer integral expressed as Uj is modified by the external 
field H through the Peierls phase factor. The original 
hopping integral tij defined on a square lattice is as- 
sumed to be t^j = t for the first NN pair, t' = — 0.12t 
for the second NN pair and t" = 0.08t for the third NN 
pair, which are selected to reproduce the Fermi surface 
topology in LSCO. 24 The average electron density per 
site is fixed at ~ 0.875. We studied the square vortex 
lattice with the unit cell size N r x N r where two vortices 
are accommodated. H = 2<p /a 2 N 2 with a the atomic 
lattice constant and </>o the unit flux. The field strength 
is denoted by H Nr {= l/N 2 ). 

The large AFM moment at underdoping decreases 
with increasing hole doping, and vanishes at overdop- 
ing due to the absence of the nesting of the Fermi sur- 
face. This magnetic tendency is simulated by changing 
U in our calculation since the following results domi- 
nantly depend on the amplitude of AFM, rather than 
the Fermi surface shape. We select U/t = 0.0, U/t = 2.6 
and U/t = 3.0, and set V/t = —2.0. These U values are 
designed to mimic the behaviors of overdoped, interme- 
diate and undcrdoped cases, respectively, in LSCO. For 
U = 0, A = l.Ot, and T c ~ 0.41i. 

To determine the eigenenergy E a and the wave func- 
tions u a (r,) and v a (ri) at the i-site, we solve the 
Bogoliubov-de Gennes equation given by 

V ( K h iJ Di ' j \ ( Ua ( r ^ \ - f ( Ua ( r ^ \ 

Y V D l - k Uj J V Mrj) ) *"* \ v a (n) ) 

(2) 

where K a ^j = —tij +Sij(Urii^ a — ji) with the chemical 
potential /U, Dij = VJ2g ^i,j^j,i+e an d a is the index of 
the eigenstate. The self-consistent condition for the pair 
potential is A hj = -|E a !( «( , '') l 'a( f, ]) tailll ( £ a/2T). 
The so-called Dopper shift effect is also included through 
phase winding of A around a vortex. The charge 
densities are n ljT = J2 a \ u a(fi)\ 2 f(E a ) and rijj = 

£>a(ri)| 2 (l-/(£7 a )). 

We calculate thermal conductivity k following the 

standard linear response theory. All the details are de- 
scribed in our previous paper studying the U — case. 22 



To recapitulate briefly, the heat current flows in response 
to a small temperature gradient. The local thermal con- 
ductivity n{r) is written as 



r' 

(3) 

where the heat-current correlation function Q xx is ex- 
pressed by the above eigenvalues and cigenfunctions. 
The spatially averaged thermal conductivity t^xx — 

J2 r K xx(i") is observed experimentally. 

In Fig. 1, we show our results of k(H)/k(Q) as a func- 
tion of H for the selected temperatures at U = (a), 
U/t = 2.6 (b) and U/t = 3.0 (c). It is seen from Fig. 1(a) 
corresponding to the overdoping case that at the lowest 
T, k(H ) increases at low H corresponding to the fact that 
ZEDOS at the core grows with H. At high T (T > T*), 
k{H) decreases with H (with increasing vortex number) 
because the heat flow n{r) is suppressed at the core. 22 
Near T <~ T* , n{H) becomes almost independent of H, 
exhibiting a plateau behavior. These behaviors are seen 
in the data x = 0.17 and x = 0.22 obtained by Sun et 
al. 18 where the corresponding T*'s are ~ 3 K and 5 K, 
respectively (see Fig. 2 in ref. 18). Various k(T) curves 
for different fields pass through 7 1 *. 25 This focal point 
temperature T* is determined by the energy width of 
the zero-energy peak at the core. This plateau behavior 
of k{H) is reminiscent of the data obtained by Krishana 
et al. 26 who attributed this to the field-induced transi- 
tion from a c4;2_ y 2-wave pairing to a pairing with a full 
gap such as d x 2_ y 2+id xy or d x 2_ y 2+is. 

The most intricate and interesting case is shown in 
Fig. 1(b) where we depict the results in U/t = 2.6 corre- 
sponding to intermediate dopings. At the lowest T, k(H) 
initially increases at low H. Then upon increasing T, 
k(H) begins to decrease monotonically, and therefore fur- 
ther increasing T, the initial suppression decreases. Thus, 
the overall curve k(H) is pushed up again at T = T mul 
(~ 0.06i when U/t = 2.6). This "reentrant" behavior 
is precisely reproduced experimentally by Sun et al. ls 
(Figs. 3(a) and 3(b) in ref. 18 for x = 0.08 and 0.10). It 
is noted, however, that the peak behaviors at low tem- 
peratures (Figs. 1(b) and 1(c)) are not noted in the ex- 
periment. In these cases, T m j n is 3-5K. We note here that 
Tmin decreases when U increases as explained later in Fig. 
3. It is understandable why there is no "reentrant" be- 
havior observed in the x = 0.14 data (Fig. 3(c) in ref. 18) 
because T m i n does not reach up to T < 7K in this exper- 
iment. In the U/t = 3.0 case shown in Fig. 1(c), where 
at a finite temperature, the local AFM order centered 
around a vortex core appears, n(H) strongly decreases 
with H . 

In Fig. 1(d), we display the field dependences of 
k(H)/k(0) for three U values at a fixed temperature. 
Depending on U, n(H) decreases or increases at low H. 
These contrasting behaviors correspond to the principal 
observation by Hawthorn et al. 19 Namely, for the un- 
derdoped samples (x = 0.06 and 0.09), dn(H)/dH < 0, 
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Fig. 1. Field dependence of k(H)/k(0). U/t = 0.0 (a), 2.6 (b), 3.0 (c) and (d) T = 0.02t for U/t = 0.0, 2.6 and 3.0. The horizontal axis 
denotes Hj^ r = l/N^ (N r the magnetic unit cell size). Note that T m ; n ~ 0.06t in (b). 



while for x = 0.17 and 0.20, dn(H)/dH > 0. This ob- 
servation is basically consistent with that by Sun et al. 18 
where the increasing and decreasing tendencies in n{H) 
is divided at the optimal doping (Fig. 4 in ref. 18). 
Hawthorn et al. interpreted this as coming from the field- 
induced thermal metal-insulator transition because the 
undcrdoped systems exhibit lack of electronic heat con- 
duction. Here, we succeed in reproducing this "transi- 
tion" by the introduction of the AFM local order in the 
undcrdoped case, which effectively removes the quasipar- 
ticles around a core, leading to poor thermal conductiv- 
ity. 




Fig. 2. Spatial structure of effective local thermal conductivity 
n(r) at H Nr = H 32 at T = O.Olt [(a) and (b)] and at T min [(c) 
and (d)], when (7 = [(a) and (c)], and U = 3.0 [(b) and (d)]. 
The vortex cores arc located at a corner and the center of the 
unit cell 



To understand this field-induced transition more 
closely, we analyze the spatially dependent thermal con- 
ductivity indicating where the heat current is 
mainly carried spatially in the vortex lattice. At low tem- 
perature, k(H) reflects the zero-energy LDOS. 22 There- 
fore, when U = [Fig. 2(a)], n(r) has a peak at the 
vortex core, where the ZEDOS peak appears. However, 
when AFM appears and ZEDOS is suppressed at the 
vortex core, n{r) is also suppressed at the core [Fig. 



2(b)]. Roughly speaking, «(r c ) at the farthest site r c 
from vortices is nearly equal to k(H = 0). Therefore, 
since the spatial average of n(r) is larger (smaller) than 
k(t" c ), we obtain k(H)/k(0) > 0(< 0) in the case of Fig. 
2(a)(Fig. 2(b)). In Figs. 2(c) and 2(d), we show n{r) at 
Tmin ~ 0.07t. In this high-temperature case, since the 
contribution to n(r) comes from the LDOS at energys 
near A , n{r) is suppressed at the core. 22 That is, the 
vortex core pushes aside the thermal flow. Therefore, the 
spatial average of k(t) is smaller than K(r c ), resulting in 
k{H) < k(0) at a high temperature. 

In the presence of AFM, the LDOS is different be- 
tween up-spin electrons (N-\ (r, E)) and down-spin elec- 
trons (N±(r, E)), depending on the sign of the AFM mo- 
ment. When N-[(r,E) > A i (r, J B) at a site, Af (r, E) < 
N^(r,E) at the nearest site. As the electron spin is 
conserved in the electron hopping, n(r) is associated 
with the smaller LDOS y /N- t (r,E)N l (r,E) , rather than 
(Ni(r,E) + Ni{r,E))/2. With increasing U, the differ- 
ence between N^{r,E) and Ni(r,E) is enhanced. This 
suppression of y/N^r, E)Ni(r, E) at E ~ A is associ- 
ated with the fact that the suppression of n(r) near the 
vortex core becomes eminent in the U/t — 3.0 case [Fig. 
2(d)] compared with the U = case [Fig. 2(c)]. 

As mentioned before, there are two characteristic tem- 
peratures T* and T m i n that govern the k(H) behavior in 
Fig. 1. T* is the temperature at which the fZ-dependence 
of k(H) is almost constant at low H and T m - m is the 
"reentrance" temperature. In Fig. 3, we show the In- 
dependence of k(H 2 4,T)/k(0,T) for the three U values 
together with the magnetic moment at the core. Since 
T* is the temperature corresponding to k(H)/k(0) = 1, 
as U increases, T* decreases and ultimately T* — > for 
U/t — 3.0. This explains why there is no plateau curve 
observed until the very low T in Fig. 1(c). 

As for the "reentrant" behavior in k(H), T m ; n occurs 
just above T*. For T* < T < T min , the slope of the 
decreasing k(H) at low H becomes more steep as T in- 
creases. However for T > T min , the slope gradually be- 
comes gentle. This corresponds to the "reentrant" behav- 
ior in k(H) as shown in Fig. 1(b). It is noted from Fig. 
3 that as T further increases, the two «;(i?24, T)/n(0, T) 
curves for U ^ merge with that for U = because 
the local moment vanishes gradually as shown in the in- 
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Fig. 3. T-dcpcndcncc of normalized thermal conductivity 
«(H 2 4,T)/k(0,T). H]v r = H 2 4(= 0.001736), [7/i=0.0(+), 
2.6(0) an d 3.0(H). The inset shows the T-dependence of the 
magnetic moment at the core site for U/t = 2.6 and U/t = 3.0. 



set of Fig. 3. It is also noted from Fig. 3 that T m ; n de- 
creases as U increases because the AFM pushes down the 
k(H2a, T)/k(0, T) curve as a whole. As mentioned above, 
this precisely corresponds to the data obtained by Sun 
et al. 18 

Let us examine the present results in a wide perspec- 
tive. As mentioned above, the experimental evidence for 
the field- induced AFM local order is abundant. Neutron 
experimental results on LSCO 1-4 are directly connected 
to the present thermal conductivity calculation. Namely, 
the observed enhancement of neutron signals under H is 
due to the local AFM order which leads in turn to the 
suppression of the zero-energy density of states around 
a core. This ZEDOS otherwise piles up as H increases. 
This suppression of ZEDOS itself has been known for 
some time in Bi2212 15,17 and YBCO systems. 16 Because 
of this suppression, k(H) decreases with H for under- 
doped cases, while n{H) increases with H for overdoped 
cases. According to Fig. 4 in ref. 18, the reversion is 
noted to occur at the optimal doping x — 0.16. This dop- 
ing dependence coincides with the observed tendency in 
LSCO by neutron scattering experiments performed ex- 
tensively for wider dopings under zero field 27 where an 
elastic peak for static AFM vanishes at overdoping. It is 
natural to expect that a magnetic field could strengthen 
this tendency. 2,3 

It should be noted that, in the mixed state, the In- 
dependence of k cannot be a simple power law T a with 
an integer a, that is, neither a = 2 for the d-w&ve zero- 
field case nor a = 1 for the high field limit because LDOS 
is spatially nonuniform (see ref. 10 for details). 

Thermal conductivity is one of the direct probes for 
observing the low-lying quasiparticles formed around a 
vortex core through its field and temperature depen- 
dences. The complex, yet rich n{H) behaviors observed 
by Sun et al. ls and Hawthorn et al. 19 yield a vivid pic- 
ture of the spatial spectral structure of low-energy quasi- 
particles through our analysis. Together with our pre- 
vious analysis 10 of site-selective NMR experiments on 



YBCO 12 and T12201 systems, 13 - 14 the present calcula- 
tion which yields a qualitatively consistent explanation 
for k(H)/k(0) unambiguously demonstrates that, (1) in 
underdoped LSCO, the local AFM order or stripe order, 
depending on the interaction strength U or doping level, 
must be present. (2) In overdopings, the ordinary (i-wave 
vortex core picture is applicable. 
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